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The role of photon scattering in shaping the lightcurves 
and spectra of 7-ray bursts 
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ABSTRACT 

We analyze the power spectra of the lightcurves of long gamma-ray bursts, dividing 
the sample in bins of luminosity, using the recently discovered variability-luminosity 
correlation. We find that the value of the variability parameter strongly correlates 
with the frequency that contains most of the power in the burst comoving frame. We 
compute the average power spectra in luminosity bins. The average power spectrum 
is well described by a broken power-low and the break frequency is a function of 
the variability parameter, while the two slopes are roughly constant. This allow us 
to conclude that scattering processes do not play a relevant role in modelling the 
lightcurves. We finally discuss in which conditions scattering may still play a relevant 
role in shaping the spectra of GRBs. 
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1 INTRODUCTION 

The study of the BATSE 7-ray bursts (GRBs) lightcurves 
has recently gained new interest thanks to the discovery of 
the variability-luminosity (Fenimore & Ramirez-Ruiz 2000; 
Reichart et al. 2001) and lag-luminosity (Norris et al. 2000) 
correlations (see also Chang, Yoon & Choi 2002). These en- 
able to assign a tentative redshift to BATSE GRBs, allowing 
to perform spectral and temporal analysis of the lightcurves 
in the burst comoving frame, where their properties are more 
closely linked to the physics of the burst itself. As an exam- 
ple, the use of these correlations enabled the possible discov- 
eries of an evolution of the luminosity function with redshift 
(Lloyd-Ronning, Fryer & Ramirez-Ruiz, 2002) and of a cor- 
relation between the peak photon frequency with the lumi- 
nosity of the GRB (Lloyd-Ronning & Ramirez-Ruiz 2002). 

The variability-luminosity correlation predicts that the 
peak luminosity of a burst is correlated to its degree of vari- 
ability, whose operational definition is related to the nor- 
malized variance, or the root mean square of the deviations 
from a smoothed version of the light curve. More variable 
lightcurves have larger luminosities. The lag-luminosity cor- 
relation is instead based on the measure of a temporal lag 
between the detection of high energy and low energy pho- 
tons. It is found that the more the burst is luminous, the 
smaller is the lag. Besides being extremely useful tools, these 
correlations also call for an explanation of their origin. Sev- 
eral possible interpretations have been discussed in the lit- 
erature. In particular, it is shown that an underlying corre- 
lation between the isotropic equivalent luminosity and the 



Lorentz factor of the fireball can explain both the corre- 
lations (Ramirez-Ruiz & Lloyd-Ronning 2002; Kobayashi, 
Ryde & MacFadyen 2002; Salmonson 2000; Meszaros et al. 
2002). In some of these works, however, a significant role is 
attributed to the modification of the temporal properties of 
the light curves as a consequence of scattering by cold or 
hot electrons (see also Panaitescu, Spada & Meszaros 1999; 
Spada, Panaitescu & Meszaros 2000). The role of scattering 
seem to be strengthened by the need of explaining the peak 
photon frequency vs. isotropic equivalent luminosity corre- 
lation (Lloyd-Ronning & Ramirez-Ruiz 2002), which is not 
naturally predicted in the internal shock scenario (Ghisellini, 
Celotti & Lazzati 1999). 

In this paper, we concentrate on the variability lumi- 
nosity correlation, using power spectra as a diagnostic to 
investigate several related issues. First (§2), we compute the 
dominant frequency of the power spectrum for each GRB in 
our sample of 220 lightcurves. We find that this dominant 
frequency strongly correlates with the variability measure. 
We then divide the sample in bins of variability and com- 
pute the average power spectrum. We find (§3) that the 
power spectra are self-similar (Beloborodov, Stern & Svens- 
son 1998; 2000; hereafter B98 and BOO; Chang & Yi 2000) 
in all the variability bins, with a break at a frequency that 
correlates with the variability parameter. No sign of a cut- 
off in the spectrum at large frequencies is observed. We (§4) 
develop a shot-noise model for the lightcurve, taking into ac- 
count the effect of scattering in Fourier space (§ 5). In §6 we 
discuss our results, showing how scattering processes cannot 
be responsible for the variability-luminosity correlation, and 
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Figure 1. The dominant frequency in the burst comoving frame 
vs. the burst variability parameter. Solid points show the bes 
class events (see Sect. 2 for details) while gray symbols show the 
lower class events. A correlation is clearly evident. The dashed 
line shows a relation v$ oc V 2 as a guideline for the first class 
events. Vertical dotted lines show the boundaries of the subclasses 
in which the sample has been divided. 



also constraining the regions of the parameter space where 
photon scattering on cold electrons may imprint detectable 
signatures on the photon spectra of GRBs without leaving 
a detectable trace in their power spectra. 



2 DATA ANALYSIS 

We have adopted the whole sample of GRBs for which 
Fenimore & Ramirez-Ruiz (2001) computed the variabil- 
ity parameter and a tentative redshift. For each lightcurve, 
we have computed the power spectrum (hereafter PSD) as 
PSD(v) = |/(f)| 2 , where f(u) is the Fourier transform of 
the lightcurve. Lightcurves were binned to 64 ms resolution 
and considered from the BATSE trigger time for a total du- 
ration of t — 3 * Tgo , where T90 is the time containing 90% of 
the total burst emission. A cubic function was fitted to the 
time interval before and after the one considered above in 
order to remove the background. The effect of background 
removal was however weak, affecting only the smallest fre- 
quencies. 

The dominant frequency Vd was detected as the fre- 
quency at which the power per unit decade frequency is max- 
imized (something similar to the peak of the v F{v) photon 
spectrum), i.e. the maximum of the function v PSD(v). We 
adopted this definition since it does not involve any fitting 
procedure and is therefore the most objective. The definition 
of the dominant frequency must confront, however, with two 
problems. First, we must consider that the time series is lim- 
ited (see below), secondly the presence of a white noise com- 
ponent, which is due to the Poisson noise of the lightcurve 
(which is inevitably affected by photon count statistics). 
This component results in a flat power spectrum (see e.g. 



Leahy et al. 1983 and references therein) which produces 
a false peak at the highest considered frequency when the 
v PSD(v) function is considered for lightcurves in which the 
signal is weaker than a certain threshold. In order to avoid 
this contamination of spurious peaks, the spectra were all vi- 
sually inspected, and every time the peak frequency was in 
the flat part of the PSD, the considered frequency range was 
shrunk in order to recompute the maximum in a shorter fre- 
quency range, where the signal-to-noise ratio is larger. This 
procedure, unavoidably, reduces the objectivity of the analy- 
sis, and for this reasons the lightcurves in which the peak se- 
lection was modified by hand were nagged. In Fig. ^ we show 
the comoving frame dominant frequency vs. the variability 
parameter for all the best class lightcurves, in which human 
intervention was either not required or clearly unbiased. A 
clear correlation is present. The correlation, according to the 
correlation coefficient statistics, is highly significant, with a 
vanishing probability of being spurious P — 3 x 10~ 25 . It 
may be argued that the correlation is artificially made by a 
correlation between the redshift and the variability (which 
is indeed present in the data). In order to test for this to 
be true, we computed the above probability also for the ob- 
served frame frequency. Again, a strong correlation is found, 
with probability of being spurious P = 8 x 10 -12 . The lower 
class bursts, in which human intervention was necessary and 
possibly biased, or in which the dominant frequency was the 
lowest one, are plotted in light gray in Fig. It is possible 
to see how these points do not invalidate the correlation, 
even though the dominant frequencies tend to be, as ex- 
pected, larger than for the best cases. In the following we 
will perform our analysis only on the subsample of 143 best 
lightcurves shown in Fig. even though extending it to the 
whole sample does not affect significantly the conclusions. It 
is worth stressing that such a correlation is not entirely sur- 
prising, since it simply says that the variability parameter 
has something to do with the mean frequency over the PSD. 
What is relevant, as we are going to discuss in much more 
detail in the following, is that this correlation tells us that 
a change in the variability parameter V reflects a change in 
the whole PSD rather than the suppression of high (or low) 
frequency power. 

When computing the PSD of an experimental time se- 
ries, one has to take into account all the effects due to the 
data binning and to the finite duration of the data. In partic- 
ular for GRBs, this second issue is particularly tricky, since 
the data are truncated not by hand, but by the intrinsic du- 
ration of the event. Consider an infinite function foo(t). The 
binning process can be modelled as the convolution with a 
square filter b(t) and the multiplication with a lattice s{t) 
(a series of Dirac 5 functions) with spacing equal to the full 
width of the square filter. The truncation of the data is rep- 
resented as the multiplication with a window function w(t) 
which is usually assumed to be the characteristic function 
of the considered interval, but may have different shapes in 
the case of GRBs (see below). Thanks to the properties of 
multiplication and convolution in Fourier space, the Fourier 
transform of / can be written as (Lazzati & Stella 1997; 
hereafter LS97): 

/>) = {[/ocM*™H] l{y)}*ii{y) (1) 

where the symbol * indicates a convolution. Let us an- 
alyze the implications of Eq. [l] from right to left. The 
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Figure 2. Average power spectra of burst in the six variability subclasses denned in the text (see also Fig. |l]). The spectra have been 
binned in frequency in order to have a constant number of points per logarithmic unit frequency. The gray line shows the best broken 
power law (plus constant) fit. Error bars are derived from the dispersion of the sample. 



Fourier transform of the lattice of spacing St (in our case 
64 ms) is the reciprocal lattice, i.e. a lattice with spacing 
2n/St ~ 100 s. This term is important if periodicities are 
present in the lightcurve, since it gives rise to the well known 
phenomenon of aliasing. In our case, since the power spectra 
are monotonically decreasing with frequency, it is irrelevant. 
The Fourier transform of the binning function has a form 
(see, e.g., LS97) 



b(v) oc 



sin (tt St v) 
7T St v 



(2) 



This is the most dangerous term, since it introduces a break 
at a frequency v ~ 1/St ^ 16 Hz in the observer frame. The 
largest of our dominant frequencies, in the observer frame, 
is maxf^) = 0.5 Hz. We conclude that this effect does not 
play a relevant role in our analysis. Finally, the effect of the 
window function is to smooth the observed PSD. In the most 
common case of a square function, the smoothing kernel is 
similar to the function in Eq. ^| but with a smaller width. 



Again, this term is important is small scale features overlaid 
on the spectrum are concerned. For any reasonable shape of 
the window function we can ignore its effect. 



We therefore conclude that the correlation of Fig |l] is 
real. It is worth noting that the typical frequencies, that 
seem to dominate the determination of the variability pa- 
rameter V , are small, even in the comoving frame. This sug- 
gests that it is the ~ 1 second variability which is physically 
linked to the luminosity of the GRB rather than its shorter 
time scale fluctuations (see also below) . There is also a sug- 
gestion, in the data, that the very small V bursts belong to 
a different family in terms of their temporal properties. We 
will discuss this in the following, finding that the average 
power spectra seem to be more consistent with a continu- 
ation of the average spectral properties rather than with a 
different subclass. 
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Table 1. Results of the fit of the average power spectra in Fig. g 
The rightmost column show the significance for the existence of 
the break. 



3 AVERAGE POWER SPECTRA 

We have divided the sample of 143 GRBs in six bins of the 
variability parameter V . Each bin contains ~ 25 lightcurves, 
with the exception of the smallest V bin, in which only 12 
lightcurves are contained. This difference is due to the need 
of keeping these bursts isolated since they seem to belong 
to a different class. The lightcurves, background subtracted, 
were normalized to contain the same numbers of photons 
and the power spectra averaged at the same comoving fre- 
quency. The errors were derived from the dispersion of the 
sample. 

This process is similar to what was done by B98 and 
BOO. Their average power spectrum presents however three 
fundamental differences. First, they performed the average 
at the same observed frequency, instead of the comoving 
frequency. This was due to the fact that redshift estimates 
were not available at that time. Secondly, they decided to 
normalize their lightcurves to the same peak photon lumi- 
nosity rather than to the same photon fiuence. We decided to 
use this second normalization since the scatter in the sample 
of spectra was lower in this case (see B98 for a discussion). 
Finally they built the average of all the lightcurves, while 
we subdivide our sample in bins of variability. 

The resulting average power spectra are shown in Fig. ^| 
The variability parameter increases from left to right and 
from top to bottom. The spectra seem to be made by a 
broken power law, sinking into the white noise at large fre- 
quencies. It is worth to note that our normalization criterion 
does not produce uniform white noise values. For this reason 
the subtraction of this component from the average spectra 
is not possible, and it is also dangerous to investigate the 
frequency region where the average spectra are dominated 
by white noise. 

We fit the average spectra with smoothly a broken 
power-law function of the form 



PSD(u) 



2F 



1/2 



+ K 



(3) 



where ai is the slope at frequencies smaller than the break 
Vq, a.2 is the slope at larger frequencies and K is a constant 
that takes into account the contribution of the white noise. 
Fq is the value of the PSD at the break frequency. 

The results of the fit are reported in Tab. jj] and in Fig. ^| 
Several interesting remarks can be derived from the results. 
First, all the power spectra are well fit by the model. More- 
over, the slopes of the broken power-law before and after the 
break are roughly consistent with being constant, with the 
break frequency vo being the only quantity that evolves with 
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Figure 3. The break frequency and high energy slopes of the 
average power spectra versus variability. The upper panel shows 
the knee frequency as a function of variability. The first point 
(open circle) is for the lowest variability class. The lower panel 
shows the high frequency spectral index versus the variability 
parameter. 



the variability parameter V. A marginal evidence of evolu- 
tion of the high frequency power law index is also present, 
especially if the smallest variability bin is included in the 
sample. A second important remark is that the high fre- 
quency slope is larger than what found by B98. This differ- 
ence is not due to the redshift correction nor to the different 
normalization. In fact, if we average among them the six 
power spectra of Fig. ^, we recover with good accuracy the 
—5/3 slope in B98. It looks therefore like that the slope can 
be attributed to the convolution of the break frequencies 
in different variability bins. Finally, even though the lowest 
variability bin values are different from the rest of the sam- 
ple, it is not possible to define it as a separate class given the 
present data. The apparent segregation of the dominant fre- 
quencies in Fig. |l| may be due to the fact that the dominant 
frequencies of the very low variability bursts approach the 
lower cut in the analyzed frequencies, and include therefore 
an additional noise term. 



4 SHOT NOISE MODEL 

In the internal shock model for GRBs (Rees & 
Meszaros 1994) the lightcurve is modelled as the random 
superposition of a number of pulses with similar shape 
stretched in shape and scaled in luminosity according to 
some prescription for the ejection of shells by the inner en- 
gine (Kobayashi Sari & Piran 1997; Panaitescu et al. 1999; 
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Spada et al. 2000, 2001; Ramirez-Ruiz & Lloyd Ronning 
2002). This lightcurve is similar to the shot noise model 
supposed to play a relevant role in the red-noise component 
observed in the power spectra of X-ray pulsars (LS97 and 
references therein). 

Consider a normalized pulse of shape p(t — to, T r ), with 
rise time r r and unit decay time and fiuence peaking a time 
to- A GRB lightcurve can be described as a random super- 
position of TV such pulses: 



t-u 



(4) 



where the pulse fiuence fj, the stretching factor ctj and the 
peak time tj are randomly selected according to some pre- 
scription (for example simulating the hierarchical shock evo- 
lution of an inhomogeneous flow). The Fourier transform of 
the above equation can be written (up to a normalization 
factor) as: 



lc[u) oc ^2 f( a ^)e 



(5) 



When we compute the average power spectrum of a sample 
of lightcurves, under the assumption that each lightcurve 
is a random realization of the same underlying process, we 
compute the ensemble average of the square of the modu- 
lus of Eq. |H| This, if the peak times of the pulses are not 
correlated with their properties and are uniform in timeQ is 
given by (LS97): 



(PSD) oc (f 2 (auj)) 



(6) 



i.e. the average power spectrum is independent on the time 
history of the pulse ejection and on the pulse fiuence distri- 
bution, but depends on the distribution of the pulse dura- 
tions. To understand the effect of the above conclusion, we 
consider a simplified exponential pulse profile: 



P(t) = e *X(o,oo) 



(7) 



where X(a,b) is the characteristic function of the interval 
(a,b). In addition, we consider a power-law distribution^] of 
a values between a minimum value a m and a maximum 
value olm'- 

n(a) oc ct a X(a m ,a M ) (8) 

The power spectrum of the single pulse of Eq. ^ is a 
Lorentzian function, i.e. (roughly speaking) a constant for 
v < 1/2-k and a power- law u~ 2 for u > 1/2tt. For the 

1 It is a well known result that the average width of the pulses 
does not evolve during the bursts (Ramirez-Ruiz & Fenimore 
2000), but their frequency and/or fiuence may be larger in the 
early phase of bursts. This correlations takes the form o a non- 
square window function in Fourier space, and are not relevant 
here, as discussed in § 2. Finally, it may be wrong to assume that 
the pulse fiuence does not depend on the pulse duration (Ramirez- 
Ruiz & Fenimore 2000). Should this assumption be wrong, it 
would reflect in a different numeric value for Eq. 

2 It has been found that the duration distribution of "well 
separated" pulses is well described by a Log-normal function 
(McBreen et al. 1994; Li & Fenimore 1999). We prefer to adopt a 
power-law distribution here since it fits better our average power- 
spectra and since our PSDs do include non well separated pulses. 
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Figure 4. Power spectra for stretched exponential pulses as de- 
fined in Eq. [121 for the set of parameters T r = 3 s, = 30 s, and 
the values of s indicated. The dotted line shows the power spec- 
trum for a single exponential profile while the dashed line shows 
the PSD for a double sided exponential profile. 



stretched pulse, the PSD is flat for v < 1/(2ttq) and a power 
law afterword. Consider now the average PSD of Eq. ^| for 
the a distribution defined in Eq. 0. If a < — 1 the PSD 
will be dominated by the shortest pulse, and will retain 
a Lorentzian shape. For a > 1, it will be dominated by 
the longest pulse, again retaining the Lorentzian shape. For 
— 1 < a < 1, a new power-law branch will appear, in the 
range 1/(27tqm) < v < l/(27ra m ), with a slope v~^ a+1 \ 
Since for a > 1 the PSD is dominated by the longest pulse, 
here and in the following we will define the "shortest rele- 
vant pulse" as the shortest pulse that have an influence in 
the PSD shape. This pulse may not be the shortest observed 
in the lightcurve. For example, consider a broken power-law 
distribution of pulse durations, with — l < a < 1 at short 
durations and a > 1 at short durations, the break in the 
power spectrum will be related to the pulse at the break of 
this distribution. We call this the shortest relevantpulse. 

The shape of the fundamental pulse in Eq. [] can be 
made more complicated. Consider as an example, a double 
exponential pulse, of functional shape: 



p(t,T r ,T d ) 



1 



Tr + T. 

Its PSD can be written as 



1 +LU^ {UJ^T? + 1 + T?/T3) 



0) 



(10) 



where lo = 2ivu is the pulsation. If r r <C Td, Eq. [10| can be 
approximated as a double broken power-law: 



|p(w)| 2 CL I V 2 



V < 1/271-Td 

l/27rr d < v < 1/27TTV (11) 

[ V > l/27TT r 

Again, the ensemble average of Eq. ^| can add different 
power- law branches if — 2 < a < 0. Finally let us consider 




Figure 5. Transfer function shapes for a uniform cloud of perfect 
scattering centers with radius R and opacity tj-. Different line 
styles show different opacities. For tt ^S> 10 the shape of the 
transfer function becomes scale invariant. Each curve is obtained 
by a Monte Carlo simulation with 10 5 input photons. 



Figure 6. Power spectral density (normalized to the value at zero 
frequency) of the transfer functions shown in Fig. tj, as a function 
of the dimcnsionless frequency (Rtt/c)ttv. For tt > 1 the PSDs 
show a prominent break at the unit value of the dimensionless 
frequency. 



a stretched double exponential (Norris et al. 1996) of equa- 
tion: 



p(t,T r ,T d ) OC 



,o) + e 



-(*)' 



X(0, + oo) 



(12) 



For s < 1 this gives a pulse with a very spiky core and broad 
wings, while for s > 1 the resulting pulse has a square shape. 
The PSD cannot be analytically computed. In Fig. ^ we show 
the shape of the PSD for a set of values of s, compared to 
the PSD of the single and double exponential pulses. For 
s < 1 the effect is of smoothing out the breaks and extend- 
ing the power-law branch u~ 2 , while the effect of s > 1 is 
more complex. For 1 < s < 2 the break frequency is moved 
to larger values, while the slope of the power-law decay is 
increased. For s > 2 the PSD starts to deviate from the 
simple form, with the appearance of "absorptions". In an 
ensemble average these small scale features will be erased 
and we are left with the envelope, which yields again a very 
steep power-law slope. 



5 PHOTON SCATTERING 

It has been proposed that photon scattering may play a 
relevant role in shaping both the temporal (Panaitescu et 
al. 1999; Spada et al. 2000; Kobayashi et al. 2002; Ramirez- 
Ruiz & Lloyd-Ronning 2002) and spectral (Meszaros & Rees 
2000, Meszaros et al. 2002) properties of GRB lightcurves. 
The effect of photon scattering has a distinctive signature 
in Fourier space, and the power spectrum of the lightcurves 
is then the best way in which the importance of scattering 
can be evaluated. 

Consider a source of photons producing a flash (ideally 
a Dirac 8 function in time) in the centre of a cloud of ideal 



scattering particles^]. Even though the central source emit- 
ted a photon impulse, an observer located outside the cloud 
would detect a light pulse with a finite duration and a par- 
ticular shape. This is due to the fact that different photons 
make different paths, finally reaching the observer after be- 
ing scattered many times in random directions. The shape 
of this pulse will be function of the opacity of the cloud to 
photon scattering only, up to a scale factor that takes into 
account the size of the cloud. 

Consider now the central source itself producing a signal 
with finite duration and with a given shape s(t). This can 
be approximated as an infinite series of spikes, each of them 
being detected by the outside observer as a pulse of shape 
k(t). In mathematical terms, the observer will detect a signal 
S(t) given by the convolution of the original function s(t) 
times the transfer function k(t): 



S(t) = s(t) * k(t) 



s(t~t')k(t')dt' 



(13) 



Thanks to the properties of Fourier transforms, the PSD of 
the detected signal 5* is the product of the spectra of s and 
k. 

In the case of photon scattering by free electrons, we 
consider a cloud of radius R and uniform density n, with 
opacity r = Rugt- We neglect the angular dependence of 
the Thompson cross section. In Fig. [| we show the transfer 
functions obtained by Monte-Carlo calculations of scattering 
for a set of optical depths. For r < 1 the transfer functions 
are Dirac 5 function with a small tail at large times, but 
for r > 1 the vast majority of photons undergo at least 

3 Ideal means in this context, that the photons are never ab- 
sorbed and that the scattering properties of the cloud do not 
depend on the photon energy. 
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one scattering and the transfer function becomes much more 
smooth. For r > 10 the opacity itself is not a parameter any 
more but a scale factor: 

K> w ^ = i k ^{t) (14) 

This behavior is reflected by the shape of the transfer 
functions in Fourier space. In Fig. ^|we show the PSD of the 
same transfer functions shown in Fig. |^. For small opacity, 
as expected, the PSD of the transfer function is almost a 
unit constant, so that S(u>) ~ S(ui). As the opacity increase, 
a pronounced break appears at a frequency 



For intermediate opacities, the PSD can be roughly de- 
scribed as a power-law for v > V\, while for large opacity 
the break takes the form of an exponential cutoff. In conclu- 
sion, the effect of photon scattering is to create a clear break 
in the power spectrum, suppressing all the frequencies larger 
than Vb- In Fig. [?] we show the result of a set of simulations 
that include scattering at various degrees. The left column 
shows simulated lightcurves with a single exponential shot 
noise model (Eq. ^) with 20 pulses. The pulse fiuence has a 
log-normal distribution, while the pulse duration Td is dis- 
tributed according to Eq. ^Jwith a — —1/2, a m = 10 s and 
a,M = 250 s. The top lightcurve has no scattering, while the 
second has r = 3 and the third has r = 10. In all cases 
R/c = 1. The last curve has been made by smoothing with 
t = 10 only the 7 shortest pulse. This may simulate more 
closely the effect of scattering in the internal shock scenario, 
where shorter pulses are produced closer to the inner engine, 
where the opacity of the flow is larger (see, e.g., Spada et 
al. 2001 and references therein). The second column shows 
the PSD of the lightcurves shown in the first column, while 
the third column shows the average PSD of 100 curves gen- 
erated according to the same pulse process. The difference 
between the scattered and unscattered PSDs is shown by 
overlaying the unscattered PSD with a light grey curve. 

As mentioned above, the effect of photon scattering in 
defining the duration of GRB pulses has been considered 
in a number of works (Panaitescu et al. 1999; Spada et al. 
2000; Kobayashi et al. 2002; Ramirez- Ruiz & Lloyd-Ronning 
2002). It is worth mentioning, however, that the way in 
which this was done is not completely accurate, especially 
when the power spectrum is concerned. In fact, in numeri- 
cal works, the duration of a pulse is computed as the sum in 
quadrature of the intrinsic duration (due to the shell width 
and the angular spreading times) plus the scattering dura- 
tion. This total duration is then used as a scale parameter for 
the pulse duration (the same as a in Eq. ^), assuming a fixed 
functional shape for the pulseQ. On the contrary, when the 
diffusion time scale is comparable to or larger then the in- 
trinsic time scale, the shape of the pulse should be changed, 
smoothing out its sharp features but leaving the long time 
scale features unaffected. When simulated lightcurves are 
observed in Fourier space, in both cases one has as expected 
that the average frequency has decreased. However, the real 

4 The approach of Kobayashi et al. (2002) was indeed different, 
in fact they assumed that all the collision between Thomson thick 
shells, typically those dominated by diffusion, do not produce any 
emission but reconvert the internal energy to outflow motion. 



scattering decreases the mean frequency by suppressing the 
high frequencies, while the above approximation rescales the 
whole spectrum at smaller frequencies. The result is that, 
while real scattering produces a prominent break in the spec- 
trum, this approach leaves the spectral shape unaffected. 



6 DISCUSSION 

How do the shot-noise theory compares to the average spec- 
tra derived in § 3? The observed spectra show a clear break, 
which we show is strongly correlated with the variability pa- 
rameter. The power-law slope changes from —2/3 before the 
break to ~ —2 after the break. 

Let us for the moment neglect photon scattering. The 
fact that the largest slope is ~ —2 suggests that an exponen- 
tial shot noise model can easily reproduce the observations, 
with a typical pulse decay time longer than the rise time 
(Norris et al. 1996). The fact that the low frequency slope is 
not flat, can be easily accommodated invoking a power-law 
distribution of pulse durations (cfr. Eq. ^|) 

n(r d ) oc r~ 1/3 (16) 

The correlation between the break frequency and the vari- 
ability (or the GRB luminosity) can therefore be interpreted 
as a correlation between the typical shortest relevant pulse 
duration and the luminosity. This does not mean that the 
lightcurve cannot have shorter pulses. In fact, if the distri- 
bution of pulse durations has a break, shorter pulses can be 
present but not contribute to the PSD. 

Let us now consider the effect of photon scattering. In 
a very simple scenario, all the pulses may go through a scat- 
tering screen with opacity r and radius R. In this case, a 
clear break should be present in the PSD at a frequency 
v = c/(ttRt). It is straightforward to show that the ob- 
served break cannot be due to scattering but must, instead, 
be attributed to the intrinsic properties of the unscattered 
pulses. In fact, the change of slope before and after the break 
is only of unity, while scattering would require either a jump 
of 2 in slope (intermediate opacity) or an exponential cut- 
off (r > 10). Alternatively, one can consider a scenario in 
which the opacity and size of the screen are different for 
different pulses. In order to be consistent with the lack of 
a pronounced break in the observations, this requires that 
the scattering time scale tR/c is smaller than r^. In this 
case, the PSD of the single pulse would be affected only in 
the power-law tail and its effect would not be observed in 
the average PSD. In order to explain the break in the aver- 
age PSD, however, one should consider a broken power-law 
of pulse durations. With such conditions one can explain 
the observed average power-spectrum. However, the role of 
scattering is not relevant in the shape of the PSD, and so it 
cannot be important in the measure of the variability param- 
eter V. It may be also envisaged a case in which tR/c 2> Td 
for all pulses, with a broken power-law distribution of the 
values of tR suited to mimic the average PSD shape. In this 
case, however, the lightcurve would be entirely dominated by 
scattering. The pulses should have the shape of the transfer 
functions k(t), in contradiction with what is found in their 
direct analysis (Norris et al. 1996). 

A more detailed modelling of an internal shock process 
should however take into account that different pulses may 
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Figure 7. Simulated lightcurves and power spectra of a shot lightcurve with various degrees of scattering (see text for more details). 
The leftmost column shows the a sample simulated lightcurve, whose PSD is shown in the central column. The rightmost column shows 
the average PSD of 100 realizations of such a lightcurve. The first row show the case with no opacity and scattering. In the second row 
all the pulses are smeared by a scattering slab with Rt/c = 3. In the third row all the pulses are smeared by a scattering slab with 
Rt/c = 10. In the last row, the properties of the slab are held the same, but only the shortest pulses are smeared while the longest 
are unperturbed. In the rightmost column, the vertical dashed line shows the position of the scattering break, the vertical dotted line 
shows the position of the pulse break and thegray line show the unscattered average spectrum of the uppermost row. All the simulations 
are made by single exponential shots (Eq. bj) and contain 20 pulses. The pulse fluence has a log- normal distribution, while the pulse 
duration is distributed according to Eq. Hwith a = —1/2, a m = 10 s and olm = 250 s. 



not be scattered with the same value of tR/c since the short- 
est pulses are likely to be produced closer to the engine, 
where the relativistic wind is more dense and opaque. To 
mimic such a case, we consider a lightcurve in which only 
the shortest pulses underwent scattering. In this case, the 
PSD should show a steepening break followed by a flatten- 



ing break when the unscattered pulse component becomes 
dominant (see lowest right panel of Fig. |7|). Again, this is 
not observed in real data (see Fig. ^). 

It may be argued that these models are too simplistic, 
since in real simulations each pulse is smeared with its own 
value of the parameter tR/c. What emerges from our anal- 
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ysis is that no clear signature of scattering is present in the 
PSD data. Is may still be possible that a more detailed nu- 
merical simulation for the evolution of the flow can include 
scattering in such a way that its effect is relevant but does 
not produce a cutoff in the PSD. This requires however a 
fine tuning of the distribution of the parameter tR/c and 
that the shape of at least part of the pulses is dominated 
by the transfer functions of Fig. ^. Such simulations, with 
a proper treatment of the scattering, are called for if the 
proposed link between scattering and variability has to be 
believed. 

Even though scattering processes are shown to be not 
important in determining the variability parameter V , it is 
still possible that they have some relevance in shaping the 
spectra of GRBs. In fact, if the scattering screen is small 
and thick, the cutoff frequency can be large, in a range in 
which the PSD is dominated by the white noise. This kind 
of scattering would not be influent for the temporal pulse 
profile, but the photon energy may be shifted by 

Ae ~ (4fcT - e) (17) 

where T is the temperature of the scattering electrons. This 
mechanism for modifying the typical energy of the photons 
in GRB spectra has been proposed and discussed, e.g., in 
Ramirez-Ruiz & Lloyd-Ronning (2002) and Meszaros et al. 
(2002). In particular, assuming that the comoving electron 
temperature is small and asking that the photon energy is 
sizably modified (Ae ~ e) one can show that the opacity 
must be 

r~-^«VF (18) 

where V is the bulk Lorentz factor of the flow and the right- 
most term holds if the initial observed energy of the photon 
was closed to m e c 2 . 

Before discussing this issue in detail, we consider that 
all the computation described above are relevant if the scat- 
tering screen is comoving with the relativistic flow. In fact, in 
order to preserve the burst variability, a scattering screen at 
rest in the frame of the host galaxy needs to have an optical 
depth much smaller than unity. For this reason, the screens 
that we consider here are comoving with the flow, they can 
be the same shell in which the radiation is produced, or the 
total contribution of previously ejected shells. 

In order for the scattering to be important in shaping 
the spectrum of GRBs, one have to assume a large opacity, so 
that the cutoff in the power spectrum should be exponential. 
Consider the bursts with largest variability (lower left panel 
of Fig. |^). The spectrum is well described by a power-law up 
to observed comoving frequencies v ~ 10 Hz. We can then 
constrain the comoving size of the scattering screen to be 

R' < 10 9 T 1/2 (19) 

In the framework of internal shock, the comoving width 
of the shell is given by A' = r/r ~ roT, where ro is the size 
of the shell at the moment of the ejection. The constraint |l^ 
implies then 

r £ 10 8 r" 1/2 (20) 

which is consistent with the internal shock picture only if 
the inner engine produces many shells, of the order of 1000 



shells in a burst lasting for 10 seconds. Without being con- 
fined in the standard internal shock model, one can envisage 
a scenario in which the energy is liberated in the fireball at 
radii r < roF 2 . In this case, it must be considered that the 
photons will be advected by the flow, which will become op- 
tically thin to radiation (due to the shell expansion) in a 
time scale smaller than the diffusion time scale of the pho- 
tons. In this case, the break frequency would be expected to 
be at v ~ cV 2 /(2ir R), which is always above the measured 
10 Hz. 



7 CONCLUSIONS 

We have computed the average power spectrum density for a 
set of GRB lightcurves, divided in 6 bins according to their 
variability properties and we have developed a shot noise 
model to be compared with them. Photon scattering is self 
consistently included in the model. We find that: 

• The variability parameter V as computed by Fenimore 
& Ramirez-Ruiz (2000) strongly correlates with the domi- 
nant frequency of the spectra, being defined as the maxi- 
mum of the v PSD(u) function. This frequency is relatively 
small, suggesting that the luminosity of GRBs is related to 
the variability properties of the lightcurve in the low fre- 
quency range (0.1 -f- 10 s). 

• The average power spectra are well described by bro- 
ken power laws with low frequency slope v~ 2 ^ and high 
frequency slope v~ 2 . The break frequency is a function of 
the variability parameter V . 

• The PSDs can be easily interpreted as due to the su- 
perposition of random similar shots (with a distribution of 
decay times and fluences) with double exponential shape 
(possibly stretched, see Norris et al. 1996). In this case the 
position of the break frequency (and therefore the variabil- 
ity parameter V) is related to the shortest relevant pulses 
in the duration distribution of pulses. 

• We show that photon scattering should imprint a de- 
tectable break or cutoff in the PSD. The lack of such a sig- 
nature makes untenable all models in which the luminosity 
variability correlation is ascribed to the smoothing of the 
shortest time scales in low-luminosity GRBs. 

• Under certain circumstances (a compact engine or devi- 
ations from the standard internal shock picture) it is possible 
to find models in which photon scattering plays a relevant 
role in shaping the spectra of GRBs consistent with the mea- 
sured PSD. 

To conclude, it must be emphasized that the results we 
have presented do depend at some level on the existence of 
a tight correlation between the burst variability and lumi- 
nosity. In particular, all the spectral analysis is made on the 
rest frame frequency, which is calculated thanks to the red- 
shift guessed from the above mentioned correlation. How- 
ever, the results still are valid and relevant even if such 
correlation should proved not to be real. In that case, the 
variability would not be related to the luminosity, but still 
we should conclude that scattering is not relevant in shap- 
ing the lightcurves of GRBs. To test this, we performed the 
same PSD analysis on the observed frequencies. The average 
spectral shapes do remain the same, even if the values of the 
brake frequencies are different and, again, there is no sign 
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of a cutoff that may indicate that scattering is relevant to 
determine the degree of variability of the lightcurves. 
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